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I.  INTRODUCTION 


Recent  advances  in  aircraft  and  radar  technology  have  given  rise  to 
changes  in  aircraft  construction  and  the  emergence  of  new  materials  used  in 
the  construction  of  aircraft.  One  of  these  changes  has  been  the  introduction 
of  dielectric  structures  for  use  both  inside  and  outside  the  aircraft,  replac¬ 
ing  the  metal  stuctures  traditionally  used.  The  impetus  for  the  use  of  these 
materials  often  stems  from  their  desirable  mechanical  properties  (strength, 
weight,  etc.),  but  the  desirability  of  aircraft  with  low  radar  cross  sections 
has  also  enhanced  their  attractiveness. 

The  use  of  these  dielectric  structures  in  modern  aircraft  has  introduced  a 
new  challenge  in  the  area  of  the  electromagnetic  properties  of  aircraft  that 
incorporate  these  new  materials.  These  properties  include  the  radar  response 
of  the  aircraft,  the  performance  of  antennas  mounted  on  them,  and  their  re¬ 
sponse  to  threats  such  as  Nuclear  Electromagnetic  Pulse  (NEP)  or  lightning. 
Whereas  the  scattering  characteristics  of  metal  structures  have  been  studied 
for  many  years  and  a  great  number  of  numerical  techniques  have  been  developed 
to  calculate  them,  this  is  not  true  of  dielectric  scatterers.  This  is  the 
result  of  two  primary  factors.  First,  in  years  past  dielectrics  were  of  mini¬ 
mal  importance  in  structures  such  as  aircraft.  Thus,  the  relative  neglect  of 
dielectrics  was  natural.  However,  the  second  factor  is  that  in  most  cases. 
Maxwell's  equations  are  much  more  difficult  to  solve  when  dielectrics  are  in¬ 
volved.  This  is  basically  because  fields  penetrate  dielectrics,  whereas  they 
are  simply  reflected  from  metals. 

A  dielectric  structure  currently  used  in  aircraft  that  is  of  particular 
interest  within  the  electromagnetic  modeling  community  is  the  thin  dielectric 
sheet.  These  structures  can  either  be  free  standing  (i.e.,  not  layered)  or 
layered  (e.g.,  conductor-backed).  Examples  of  the  former  are  aerodynamic 
control  surfaces  such  as  fins.  An  example  of  the  latter  would  be  a  fuselage 
constructed  using  conductor-backed,  thin,  dielectric  panels.  In  either  case, 
the  electromagnetic  modeling  of  these  structures  is  more  difficult  than  the 
corresponding  metal  structures  since  dielectrics  are  volume  scatterers,  as 
opposed  to  metal  structures,  which  are  surface  scatterers. 

This  report  studies  the  application  of  a  popular  numerical  technique,  the 
Finite-Difference  Time-Domain  (FDTD)  technique  to  the  problem  of  scattering  by 
conductor-backed  dielectric  plates,  with  cracks  in  the  dielectric.  Although 
previous  efforts  have  only  dealt  successfully  with  frequencies  that  are 
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relatively  high,  the  technique  presented  in  this  report  is  applicable  at  low 
frequencies  (i.e.,  where  the  crack  is  small  with  respect  to  wavelength).  The 
formalism  presented  requires  a  relatively  minor  change  in  the  basic  FDTD 
technique. 

The  section  to  follow  presents  a  brief  overview  of  some  electromagnetic 
numerical  techniques  as  applied  to  dielectric  structures.  This  is  followed  by 
an  introductory  section  on  the  basics  of  the  FDTD  technique.  Next  the  FDTD 
technique  is  extended  so  that  it  can  model  thin  dielectric  sheets.  This  is 
followed  by  a  discussion  of  the  extension  of  the  FDTD  technique  to  conductor- 
backed  dielectric  sheets.  Finally,  the  complete  model  of  conductor-backed 
dielectric  sheets  with  cracks  is  presented. 
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II.  GENERAL  COMMENTS  OH  NUMERICAL  TECHNIQUES  FOR  DIELECTRIC  SCATTERERS 

The  numerical  solution  of  Maxwell's  equations  has,  of  course,  been  a  topic 
of  great  importance  for  many  decades.  As  a  result,  there  are  a  great  number 
of  techniques  available  for  solving  a  wide  range  of  problems  of  interest  to 
the  electromagnetic  community.  Included  in  the  list  of  popular  numerical 
techniques  are  integral  equation  methods,  eigenvector  methods,  finite  element 
methods,  partial  differential  equation  methods,  asymptotic  techniques,  and 
radiative  transfer  methods,  just  to  name  a  few.  In  spite  of  the  seemingly 
large  number  of  techniques  available,  it  is  often  not  particularly  clear  which 
method  (if  any)  is  well  suited  to  a  particular  application.  Certainly  there 
is  no  one  technique  that  is  best  suited  to  all  situations. 

As  has  been  stated  earlier,  geometries  containing  dielectrics  are  parti¬ 
cularly  difficult  to  model  since  they  are  by  nature  volume,  rather  than  sur¬ 
face,  scatterers.  This  means  that  the  unknown  being  searched  for  (usually  a 
field  or  current  distribution)  exists  not  just  on  the  surface  of  the  scat- 
terer,  but  also  within.  This  obviously  adds  another  dimension  to  the  problem 
and  thus  increases  both  the  computational  and  storage  requirements  of  most 
numerical  techniques. 

The  presence  of  slots  or  cracks  in  a  dielectric  scatterer  poses  a  parti¬ 
cularly  difficult  challenge  since  it  combines  two  different  problems:  the 
volume  scattering  of  the  dielectric  and  the  nonlinear  field  behavior  near 
sharp  discontinuities.  In  spite  of  the  importance  of  this  particular  geome¬ 
trical  feature  within  electromagnetic  engineering,  little  progress  has  been 
made  in  the  analytical  solution  of  such  problems.  One  notable  exception  has 
been  at  high  frequencies,  where  the  Geometrical  Theory  of  Diffraction*  (GTD) 
can  be  used. 


Two  techniques  that  have  enjoyed  a  great  deal  of  popularity  within  the 
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electromagnetic  community  are  the  Method  of  Moments  (MoM)  and  the  Finite- 
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Difference  Time-Domain  (FDTD)  technique.  Both  techniques  have  proved  very 
successful  in  a  wide  range  of  problems,  and  of  course,  both  have  their  draw¬ 
backs.  Basically,  the  MoM  is  an  integral  equation  technique,  usually  thought 
of  as  operating  in  the  frequency  domain.  With  this  technique,  the  unknown 
(generally  a  surface  current  or  an  interior  field)  is  expanded  in  a  finite 
dimensional  basis  and  then  a  matrix  equation  is  solved  that  minimizes  the 
error  (generally  the  boundary  conditions)  in  the  solution. 


In  the  FDTD,  Maxwell's  equations  are  generally  solved  by  direct  integra¬ 
tion  in  the  time  domain.  Here,  both  time  and  space  are  divided  into  a  numeri¬ 
cal  grid.  The  electric  and  magnetic  fields  are  then  advanced  temporally  by 
calculating  the  line  integral  of  the  previously  calculated  fields  in  each  cell 
so  as  to  obtain  the  time  rate  of  change  of  each  field  component  in  each  cell. 
In  this  way,  time-marching  solutions  are  obtained  without  the  need  to  invert  a 
matrix. 

Although  both  of  these  techniques  have  been  used  with  success  to  solve 
problems  involving  dielectrics,  it  has  recently  been  shown ^  that  the  FDTD  has 
certain  advantages  over  the  MoM  when  dielectrics  are  involved.  These  advan¬ 
tages  stem  from  an  overall  savings  in  the  number  of  computer  operations  tuat 
need  to  be  performed  in  order  to  solve  problems  of  a  given  electrical  size. 

The  technique  has  been  continually  advanced  in  recent  years  so  that  a  wide 
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range  of  geometries  can  be  modeled.  ♦  *  ’  It  will  be  the  subject  of  the 
subsequent  sections  to  show  how  this  technique  can  be  augmented  so  that  it  can 
be  used  on  thin  dielectric  structures  with  cracks. 


4 


III.  REVIEW  OF  FDTD 


The  FDTD  technique  was  first  introduced  for  electromagnetic  problems  by 
3 

Yee.  To  a  large  extent,  this  technique  has  been  relatively  unchanged  since 
its  introduction.  The  development  in  this  section  will  be  specialized  to  the 
case  of  simple  conducting  scatterers  suspended  in  free  space.  This  discussion 
will  act  as  an  introduction  to  those  later  sections  that  consider  augmenta¬ 
tions  of  this  basic  technique  when  dealing  with  scatterers  containing  thin 
dielectric  sheets  and  the  possibility  of  cracks. 

Although  Yee's  original  formulation  of  the  technique  was  in  terms  of  the 

differential  form  of  Maxwell's  equations  (and  hence  the  name  FDTD),  it  has 
8 

been  shown  that,  in  many  instances,  a  superior  view  of  the  technique  is 
obtained  by  starting  with  the  integral  form  of  Maxwell's  equations  and  imple¬ 
menting  them  over  small  volumes  of  the  solution  space.  Since  it  is  this  inte¬ 
gral  formulation  that  offers  the  most  insight  to  the  dielectric  structures  to 
be  considered  in  the  following  sections  of  this  report,  the  discussion  to 
follow  will  be  from  this  perspective. 

Maxwell's  equations  can  be  stated  in  integral  form  as: 


4  H  •  dl  -  I  +  ^  /  eE  •  ds 
<£  E  •  dJi.  *  ~  /  p  H  •  ds 


[la] 

[lb] 


For  ease  of  reference,  the  first  equation  will  often  be  referred  to  as 
Ampere's  law,  and  the  second,  Faraday's  law.  As  can  easily  be  seen  from  these 
equations,  a  curling  magnetic  field  causes  a  time-varying  electric  field,  and 
vice  versa.  In  fact,  the  development  in  time  of  an  electromagnetic  field  can 
be  thought  of  as  the  back  and  forth  action  and  reaction  between  electric  and 
magnetic  fields.  The  FDTD  exploits  this  very  view  of  electromagnetic  fields. 

The  FDTD  approach  proceeds  to  solve  for  the  electric  and  magnetic  fields 
resulting  from  a  particular  incident  field  (or  source)  interacting  with  a 
scatterer  by  dividing  both  time  and  space  into  a  numerical  grid.  As  shown  in 
Figure  1,  the  space  in  the  vicinity  of  the  scatterer  is  divided  into  an  array 
of  rectangular  cells.  Within  each  cell,  six  components  of  the  electric  and 
magnetic  fields  are  evaluated  at  distinct  points  in  time.  Figure  2  shows  the 
points  at  which  electric  and  magnetic  fields  are  evaluated  .La  each  cell.  It 
should  be  noticed  that  electric  and  magnetic  fields  are  evaluated  at 


5 


interlaced  points  in  space.  This  interlacing  has  been  deliberately  engineered 
so  as  to  readily  calculate  the  line  integral  of  the  electric  field  about  each 
magnetic  field  point,  and  vice  versa.  These  fields  are  also  evaluated  at 
interlaced  points  in  time. 


FIGURE  1.  An  arbitrary  scatterer  immersed  in  a  rectangular  FDTD 
numerical  grid. 


y 


t 


FIGURE  2.  An  FDTD  unit  cell  showing  the  electric  and  magnetic  field 
evaluation  points 


The  FDTD  technique  begins  the  solution  of  the  fields  with  the  assumption 
that  the  system  is  starting  from  some  known  condition — most  often  an  assump¬ 
tion  of  zero  initial  fields.  Thus,  it  is  assumed  that  the  electric  fields  are 
known  at  each  lattice  point  at  t  =  0,  and  the  magnetic  fields  are  also  known 
at  each  lattice  point  at  t  "  — Given  this  assumption,  Faraday's  law  can  be 
invoked  along  the  three  perpendicular  paths  in  each  cell  to  arrive  at 
expressions  of  the  form: 

jr  E  •  dT  =  L  P(-H  )dy  dz  •  [2] 

^  I  O  t  o  ^  x 

Here,  the  path  of  integration,  C^,  is  counter  clockwise  around  the  perimeter 
of  any  cell  (in  this  case  in  the  x  =  0  plane),  and  the  surface  Sj  is  that 
which  is  enclosed  by  Cj.  This  contour  is  depicted  in  Figure  3  (which  shows  a 
y-z  plane  cut  of  the  three-dimensional  lattice).  In  this  figure,  the  electric 
and  magnetic  field  evaluation  points  pertinent  to  Equation  [2]  are  indica¬ 
ted.  Also,  Ay  and  Az  are  the  cell  dimensions  in  the  £  and  2  directions, 
respectively. 


2 


contous  used  for  advancing  the  electric  and  magnetic 
fields  in  homogeneous  regions. 


If  the  cell  In  question  contains  only  homogeneous  media,  the  evaluation  of 
Equation  [2]  proceeds  by  assuming  that  the  dimensions  of  the  cells  are  small 
enough  (generally  \/5  or  smaller)  in  each  direction  so  that  the  fields  within 
each  cell  can  be  considered  to  be  linear  functions  of  position.  Considering 


first  the  left  side  of  Equation  [2],  it  is  seen  that  it  consists  of  four  sep¬ 
arate  line  integrals.  Restricting  attention  to  the  integral  of  Ez  over  the 
left-hand  wall,  we  first  assume  that  along  this  wall  Ez  varies  as 


Ez  -  EZ(I , J,K)  +  A  •  z  ,  [3] 

where  z  is  measured  with  respect  to  the  center  of  the  cell  (i.e.,  the  Hx 
evaluation  point),  and  "A"  is  a  constant.  EZ(I,J,K)  is  the  value  of  Ez  at  the 
midpoint  of  the  left  wall  (the  indices  I,J,K  indicate  that  this  field  point  is 
associated  with  the  I,J,K'th  cell).  Thus,  the  contribution  to  the  line 
integral  of  Ez  along  the  left  wall  can  be  written  as 

/  E  dz  -  EZ(1, J,K)  Az  •  [4] 

z 

left 

wall 

A  similar  analysis  can  be  performed  on  the  other  three  legs  of  the  line 
integral  in  the  left  side  of  Equation  [2].  For  instance,  the  integral  along 
the  right-hand  path  yields  a  result  similar  to  Equation  [4],  except  that  the 
direction  of  integration  is  in  the  -z  direction  and  the  center  of  this  path 
coincides  with  the  sample  point  EZ(I,J+1,K).  Also,  the  same  type  of  assump¬ 
tion  as  shown  in  Equation  [3]  can  be  made  for  the  ^-directed  electric  fields 
along  the  upper  and  lower  legs  of  the  contour.  Thus,  the  complete  contour 
integral  can  be  written  as: 

<L  E  •  dl  =  [EZn(I,J,K)  -  EZn(I,J+l,K)]Az  [5] 

+  [EYn(I,J,K+l)  -  EYn(I,J,K)]Ay 

In  this  equation,  the  superscript  n  indicates  that  this  line  integral  has  been 
evaluated  at  time  t  =  nAt. 

Turning  attention  to  the  right  side  of  Equation  [2],  it  is  once  again 
assumed  that  Hx  is  a  linear  function  of  position  throughout  the  cell.  Thus, 
Hx  can  be  written  as 

Hx  -  HX(I, J,K)  +  A  •  z  +  B  •  y  •  '  [6] 

Here  the  positions  of  y  and  z  are  measured  with  respect  to  the  center  of  the 
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cell,  and  A  and  B  are  constants.  Since  this  cell  contains  only  homogeneous 
media,  p  is  constant  throughout  the  cell  and  Equation  [6]  can  be  substituted 
into  the  rigiit-hand  side  of  Equation  [2]  to  yield,  after  integration: 


a£  Jgj  “  <-V  ^  dz  -5*  HXd.J.K)  Ay  Az  |t  .  ^  ,  [71 

where  I  .  indicates  that  the  time  derivative  is  to  be  performed  at  time 

1 1  *  nAt 

t  *  nAt. 

The  final  step  in  evaluating  the  H  field  at  the  center  of  this  cell  is  to 
evaluate  the  time  derivative  in  the  above  equation  by  standard  central  differ¬ 
encing.  Thus: 

.  n+  V2  n-  V2 

p  HX(I, J,K)Ay  Az  lt  =  ^  ^  (I.J.K)  “  HX  (I, J ,K) ] AyAz ,  [8] 

where,  as  in  the  electric  fields,  the  superscripts  on  HX  indicate  the  times 
they  are  to  be  evaluated.  Finally,  substituting  Equation  [5]  and  [8]  into 
Equation  [2]  and  solving  for  HXn+  ^  (I^K)  yields 

HXn+  1,2  (1,J,K)  -  HXn~  1/2  (I,J,K) 

[EZn(I,J,K)  -  EZn(I,J+l,K)]  [9] 

+  Az  [EY“<I*J»K+1>  “  EYn(I,J,K)]  } 


This  is  the  final  FDTD  expression  for  evaluating  the  ^-directed  magnetic 
fields  in  each  cell.  Similar  expressions  can  be  developed  for  the  y  and  z 
components. 

Once  the  values  of  the  magnetic  fields  at  time  t  =  (n+  V2  )At  have  been 
calculated  in  each  cell,  an  application  of  Ampere's  law  in  each  cell  will 
yield  the  electric  fields  at  time  t  ■  (n+l)At.  The  procedure  for  evaluating 
the  integrals  follows  the  same  general  scheme,  with  the  exception  that  the 
contour  integrals  are  about  H  field  points,  rather  than  E  field  points,  as  can 
be  seen  in  Figure  3  (see  contour  C2).  A  representative  expression  for 
advancing  an  ^-directed  electric  field  in  a  cell  containing  only  homogeneous 
material  and  no  currents  is: 


9 


(I.J.K)] 


[10] 


EXn+1(I, J,K)  -  EXn(I,J,K) 

-  F-  t  [HZn+1/2(I,J-l,K>  -  HZn+1/2 

+  jj-  [HYn+1/2  CI,J,K)  -  HYn+l/2(I,J,K-l)]  } 

The  general  procedure  of  the  FDTD  technique  is  to  sequentially  evaluate 
the  magnetic  fields  over  the  entire  solution  space,  then  the  electric  fields 
at  the  next  half-time  step.  This  process  is  repeated  on  and  on.  In  this  way, 
time-marching  solutions  are  obtained. 

The  sources  of  the  incident  fields  in  the  solution  space  can  be  specified 
in  a  variety  of  ways.  If  the  source  can  be  specified  as  a  known  current  dis¬ 
tribution  relatively  close  to  the  source,  the  time  history  of  these  currents 
can  be  forced  directly  in  Ampere's  law  (Equation  [la]).  For  sources  distant 
from  the  scatterer,  it  is  generally  best  to  divide  the  total  fields  into  inci¬ 
dent  and  scatterered  components.  Since  the  incident  fields  are  assumed  to  be 
known  over  all  of  space  for  all  time,  the  FDTD  technique  only  needs  to  advance 
the  scattered  fields,  noting  that  in  cells  containing  free  space  they  satisfy 
exactly  the  same  Maxwell's  equations  as  do  the  incident  fields.  It  is  also 
possible  to  have  regions  in  the  FDTD  solution  space  where  total  fields  are 
evaluated  and  stored,  and  other  regions  where  only  scattered  fields  are  cal- 

Q 

culated.  This  scheme  is  popular  when  dielectric  scatterers  are  present. 

In  cells  containing  perfectly  conducting  material,  boundary  conditions 
must  be  enforced.  For  those  regions  in  which  total  fields  are  being  calcu¬ 
lated,  the  tangential  electric  fields  at  the  conductor  surfaces  must  be  main¬ 
tained  at  zero.  For  those  regions  in  which  scattered  fields  are  being 
calculated,  the  appropriate  condition  is  that  the  tangential  components  of  the 
scattererd  fields  must  be  maintained  to  equal  the  negative  of  the  incident 
tangential  fields.  For  this  latter  case,  it  can  easily  be  seen  that  the 
scattered  fields  literally  grow  out  of  the  scatterer  since  the  scattered 
fields  will  not  start  until  the  incident  field  "hits"  the  scatterer.  Only 
then  will  the  calculated  fields  begin. 

A  final  point  to  be  mentioned  here  is  the  treatment  of  the  fields  along 
the  edges  of  the  solution  space.  The  special  problem  here  is  that  in  order  to 
continue  to  update  the  tangential  electric  fields  along  this  boundary  using 
the  usual  Ampere's  law  procedure,  the  values  of  magnetic  fields  beyond  the 
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boundary  would  be  required.  If  one  proceeds  blindly  by  ignoring  this  lack  of 
information  and  tacitly  assumes  that  these  needed  field  components  are  zero, 
this  will  give  rise  to  non-physical  reflections  and  the  numerial  solution  will 
be  perturbed.  To  address  this  problem,  several  schemes  have  been  developed 
whereby  these  needed  fields  are  estimated  from  known  values  within  the  solu¬ 
tion  space.  These  techniques  yield  relatively  small  reflections,  although 
work  is  being  done  to  find  even  better  estimators.  The  technique  used  by  most 
current  investigatiors  (including  this  principal  investigator)  is  that  devel- 

Q 

oped  by  Mur  (usually  the  2nd  order  solution). 


IV.  THIN  DIELECTRIC  SHEETS 


The  analysis  of  scatterers  containing  solid  dielectric  components  is  very 
straightforward.  When  dielectric  material  is  present  throughout  cells  in 
which  total  (as  opposed  to  scattered)  fields  are  being  calculated,  it  is 
necessary  only  to  change  the  constituative  parameters  (i.e.,  n,  e,  and  a) 
within  those  cells. 

In  situations  where  a  scatterer  contains  thin  dielectric  sheets  one  can, 
of  course,  use  cell  dimensions  that  are  small  enough  to  resolve  the  sheet  as  a 
collection  of  solid  cells  of  dielectric.  This,  however,  would  generally 
require  cell  dimensions  throughout  the  problem  space  at  least  as  small  as  the 
dielectric  sheet  thickness,  thus  severely  encumbering  the  technique  both  in 
computational  time  and  storage  when  the  sheet  is  very  thin.  To  circumvent 
this  undesirable  situation,  a  scheme  whereby  the  cell  size  can  be  allowed  to 
be  much  larger  than  the  slab  thickness  is  desirable  and  will  next  be  developed 
in  the  paragraphs  to  follow. 

Figure  4  shows  a  two-dimensional  cut  of  a  FDTD  spatial  grid  containing  a 
thin  dielectric  sheet  scatterer.  The  dielectric,  of  thickness  d  and  permiti- 
vity  e,  is  represented  by  the  shaded  region.  The  electric  and  magnetic  field 
sample  points  pertinent  to  this  discussion  are  labeled  in  this  figure. 


♦  E* 

Hd 


FIGURE  4.  The  orientation  of  a  thin  dielectric  slab  within  the  FDTD 
grid,  showing  the  integration  contours  used  to  advance 
the  electric  and  Magnetic  fields  in  the  vicinity  of  the 
slab. 
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Viewing  these  field  component  locations  in  light  of  the  field-advancing  equa¬ 
tions  discussed  in  the  previous  section  (Equations  [la]  and  [lb]),  it  can  be 
seen  that  the  only  field  components  whose  calculation  will  be  directly  affect¬ 
ed  by  the  presence  of  the  slab  are  the  tangential  electric  fields  evaluated 
within  the  slab  and  the  tangential  magnetic  field  components  just  above  and 
below  the  slab. 

In  order  to  use  Ampere's  and  Faraday's  laws  (Equations  [la]  and  [lb],  re¬ 
spectively)  to  calculate  the  tangential  fields  in  the  vicinity  of  the  dielec¬ 
tric,  some  assumptions  must  be  made  concerning  the  behavior  of  these  fields. 
To  aid  in  this,  observations  were  made  of  the  fields  in  the  vicinity  of  an 
electrically  thin,  infinite  dielectric  slab,  subjected  to  a  pulsed  plane  wave 
with  various  angles  of  incidence.  It  was  found  that  the  magnetic  fields,  al¬ 
though  affected  by  the  presence  of  the  dielectric,  maintain  relatively  linear 
behavior  through  the  dielectric.  As  for  the  tangential  electric  fields,  it  is 
found  that  they  too  maintain  relatively  linear  behaviors,  with  a  small  change 
in  slope  at  the  air/dielectric  interface.  The  normal  electric  fields,  as 
would  be  expected,  exhibit  a  discontinuity  at  the  air/dielectric  interface  of 
value  equal  to  the  ratio  of  the  air  and  dielectric  permitivities. 

Given  the  above  observations,  it  is  reasonable  to  assume  that,  for  thin 
slabs  (d  <<  \),  all  magnetic  fields  and  the  tangential  electric  fields  can  be 
modeled  as  linear  functions  throughout  those  cells  containing  the  slabs.  On 
the  other  hand,  the  normal  electric  fields,  which  exhibit  a  definite  disconti¬ 
nuity,  must  be  modeled  with  this  discontinuity  built  into  the  model  in  order 
for  applications  of  Ampere's  and  Faraday's  laws  to  provide  meaningful  results. 

Turning  attention  to  the  evaluation  of  the  tangential  electric  fields 
within  the  slab,  it  is  noticed  that  the  evaluation  of  the  left  side  of 
Ampere's  law  around  a  contour  such  as  C3  will  proceed  just  as  in  the 
"standard"  FDTD  technique  since  the  magnetic  fields  are  still  assumed  to  be 
linear  functions  of  position.  The  evaluation  of  the  right  side  of  Equation 
[la]  must,  however,  be  modified  since,  although  the  electric  field  is  assumed 
to  be  a  linear  function  of  position  throughout  the  cell  volume,  e  is  not. 
Taking  this  into  account  and  integrating,  one  obtains: 
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EXn+1(I, J,K)  -  EXn(I,J,K) 


-  {  [HZn+1/2  (I,J-1,K)  -  HZn+  1/2  (I,J,K)] 

2  [1 

+  Az  [ 8,11+1/2  "  HYn+1/2  (I,J,K-1)]  }  . 


It  should  be  noted  that  the  above  formula  becomes  the  same  as  Equation  [10]  In 
the  limit  as  d+0  or  e-*^,  as  would  be  expected.  Also,  as  d  becomes  larger 
(but  still  small  with  respect  to  wavelength)  or  e  becomes  larger,  the  pre¬ 
dicted  value  of  EX(I,J,K)  becomes  smaller.  Again  this  is  to  be  expected  since 
the  natural  tendency  of  the  slab  to  lower  the  electric  field  within  it  will 
become  more  dominant  as  the  slab  thickness  or  permitlvity  increases.  A  simi¬ 
lar  equation  can  be  developed  for  Ey. 

In  order  to  advance  the  magnetic  fields  above  and  below  the  slab,  Fara¬ 
day's  law  (Equation  [lb])  must  be  evaluated  around  a  contour  such  as  C^. 
There,  the  evaluation  of  the  right  side  of  Equation  [lb]  proceeds  exactly  as 
if  the  dielectric  were  not  there  since  it  has  been  assumed  that  the  magnetic 
field  has  distributed  itself  linearly  throughout  these  cells,  just  as  in  a 
free-space  cell.  On  the  other  hand,  the  evaluation  of  the  line  integral  of 
E  around  the  perimeter  of  the  cell  must  be  handled  differently  due  to  the 
discontinuity  of  the  normal  component  of  E.  Figure  5  depicts  the  behavior  of 
Ez  as  a  function  of  z  in  a  cell  that  contains  the  dielectric  slab.  The  func¬ 
tion  that  describes  this  behavior  is 

j  ^  [EZ(I, J,K)  +  A  •  z  ]  ^|  <  z  <  ^—  + 

EZ  -  j  [12] 

l  EZ(I, J,K)  +  A  •  z  +  f  <  2  <  ^§ 


where  z  is  measured  with  respect  to  the  center  of  contour  C^. 
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Ez(Z) 


FIGURE  5.  The  assumed  behavior  along  the  t  direction  of  the  normal 
coaponent  of  the  electric  field  within  the  contour  of 
FIGURE  4. 

When  integrating  the  electric  field  described  by  Equation  [12]  in  the 
left-hand  side  of  Faraday’s  law  (Equation  [lb]),  it  is  seen  that  there  exists 
an  ambiguity  since,  although  the  value  of  EZ(I,J,K)  is  assumed  to  be  known  for 
previous  time  steps  and  obtainable  for  future  time  steps,  the  slope,  "A",  is 
not.  However,  "A"  can  be  estimated  from  the  known  values  of  Ez  on  both  sides 
of  the  slab.  Thus: 

AU,J)  -  [  EZn  (I ,  J ,  K)  -  EZn(I ,  J,  K-l  )  ]  ,  [13] 

where  A  (I,J)  denotes  the  slope  of  EZ  along  the  lattice  lines  x  *  I*Ax  and 
and  Y  ■  J*Ay.  Substituting  Equations  [12]  and  [13]  into  Faraday's  law  and 
integrating  around  the  contour  yields: 
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HXn+  1/2  (I,J,K)  -  HXn"  1/2  (I,J,K) 

+  a^t  {  l>  +  f  (-f _1)^  [EZ’,(I>J>K) "  Ezna.«j+i.K)] 

[14] 

+  [ A(I , J)  -  ACI.J+I)]  |  (-^  -1)  (|  -  Az) 

+  Ay  [EY^I.J.K+l)  -  EYn(I,J,K)} 


A  similar  equation  can  be  developed  for  H^.  If  the  slabs  become  thicker  than 
X/5,  it  is  best  to  model  them  as  solid  cubes. 

Equations  [11]  and  [14]  constitute  the  only  modifications  needed  to  allow 
the  FDTD  technique  to  model  thin  dielectric  structures. 
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V.  CONDUCTOR-BACKED  DIELECTRIC  SHEETS 


The  modeling  of  a  thin  dielectric  slab  backed  by  a  vanishingly  thin  con¬ 
ducting  sheet  is,  in  some  ways,  easier  than  that  of  a  dielectric  slab.  This 
is  because  the  tangential  electric  and  the  normal  magnetic  fields  at  the  di¬ 
electric/conductor  interface  are  known  (i.e.,  they  are  zero). 

An  FDTD  model  of  a  thin,  conductor-backed  dielectric  sheet  is  shown  in 
Figure  6.  The  dielectric  sheet  is  shown  as  the  shaded  region,  and  the  con¬ 
ducting  sheet  is  shown  as  the  solid  line  beneath  it.  A  typical  contour  used 
to  evaluate  the  tangential  electric  fields  at  the  air/dielectric  interface  is 
shown  as  Cg.  Here,  the  contour  is  shortened  in  the  z  dimension  as  compared  to 
"normal"  contours  due  to  the  presence  of  the  conducting  sheet  at  the  bottom  of 
the  contour.  Ampere’s  law  can  be  utilized  along  this  contour  to  evaluate  Ex 
along  the  air /dielectric  interface  if  assumptions  are  made  about  the  behavior 
of  Ex  and  H2  within  and  along  the  contour,  as  well  as  the  value  of  along 
the  conductor  surface. 


FIGURE  6.  The  orientation  of  a  thin  conductor-backed  dielectric 
sheet  within  the  FDTD  grid,  showing  the  integration 
contours  used  to  evaluate  tangential  electric  fields 
above  the  slab  and  magnetic  fields  below  the  slab. 


Since  both  Ex  and  Hz  have  zero  values  at  the  conductor  surface  and  will 
vary  linearly  near  and  within  the  dielectric,  Ex  and  Hz  can  be  modeled  as: 
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H 

z 


0 

HZ  (I,J,K)  [l  +  j  ] 


-Az 

2 


<  z  <  -d 


-d  <  z 


[15] 


and 


j  0  —■  <  z  <  -d 

E  «  { 

x  1 

(  EX(I , J,K)  [l  +  f  ]  -d  <  z  <  ~ 


[16] 


where  HZ(I,J,K)  and  EX(I,J,K)  represent  the  values  of  these  fields  at  the 
air/dielectric  interface.  Also,  z  is  measured  with  respect  to  the  center  of 
C^.  Equations  [15]  and  [16]  are  depicted  in  Figure  7. 


FIGURE  7 .  The  assumed  variations  along  the  f  direction  of  the 
tangential  electric  and  normal  magnetic  fields  within  the 
contour  C5  of  FIGURE  6. 


Substituting  Equations [ 15]  and  [16]  into  Ampere's  law  (Equation  [la])  and 
evaluating  about  contour  C5  yields: 
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EXn+1(I,J,K)  =  EXn(I, J,K) 


At 


[(£  "V  I  +  2d  ^  +  d)2]Ay 


[17] 


{  t  25-  +  d)2[H2n+1/2  (I,J-1,K)  -  HZn+1/2  (I,J,K)] 


Ay[HYn+1/2  (I,J,K)  -  h"*  1/2  ]  } 

In  the  above  equation,  HyC  is  the  £  directed  magnetic  field  along  the 
bottom  of  the  contour.  Note  that  this  location  does  not  coincide  with  any 
normal  H  field  evaluation  points.  This  has  occurred  because  the  conducting 
sheet  does  not  lie  along  a  grid  line.  To  evaluate  this  field,  use  will  be 
made  of  the  y  component  of  the  point  form  of  Faraday's  law: 


dE 

dE 

dH 

_ X 

z 

-u 

dz 

dx 

^  dt 

[18] 


In  order  to  evaluate  this  expression  at  the  conductor  surface,  the  deriva- 

3E 

x 

tives  on  the  left  side  of  the  expression  must  be  estimated.  Of  these,  -y  is 

the  easiest  to  evaluate  since  the  form  of  E_  has  already  been  specified  in 

dE 

x 

Equation  [16].  Thus,  r —  is  given  by 

0  z 

dE 

aT  -  i^xd.j.K)  •  [19] 

dE 

In  evaluating  the  term  jy,  It  Is  noticed  that  the  form  of  Ez  has  not 
(and  cannot)  be  as  easily  specified  as  was  Ex  since  an  electric  field  normal 
to  a  conducting  surface  need  not  have  zero  value  at  that  surface.  It  must, 
however,  approach  that  surface  with  zero  slope.  Thus,  it  will  be  assumed  that 

Az 

the  values  of  E  at  a  distance  of  d  +  y  above  the  ground  plane  are  sufficient 
dE 

2 

to  predict  ■y  at  the  conductor  surface.  Thus,  we  have: 


19 


dE  e 

ST  m  [EZ(I+1 , J,K)-EZ(I, J,K) ]  •  [20] 

£ 

Note  that  the  — —  multiplier  in  the  above  expression  arises  from  the  known 

discontinuity  in  Ez  across  the  air/dielectric  interface.  Equations  [19]  and 
[20]  can  now  be  substituted  into  Equation  [18]  to  yield 

H"c1/2  *  Uyl1/2~f'  1  I  [EZn(I+l,J,K)  -  EZn(I,J,K)]  •  [21] 

This  is  the  value  of  H^c  as  intended  for  use  in  Equation  [17]  to  evaluate  the 

tangential  electric  fields  at  the  air/dielectric  interface.  Note  that  in 

order  to  evaluate  Equation  [21],  the  value  of  HyC  at  the  previous  time  step 

must  be  stored.  This  means  that  computer  storage  must  be  reserved  for  these 

tangenial  magnetic  fields  all  along  the  conducting  sheet. 

The  evaluation  of  Hz  at  the  air/dielectric  interface  needs  no  special 

attention  since  the  integration  takes  place  for  a  fixed  value  of  z  and  the 

fields  H  and  Ev  and  E„  are  continuous  across  this  interface.  Thus,  since  the 

4  *  y 

variables  EX  and  EY  represent  the  values  of  the  field  at  the  interface,  the 
usual  calculation  of  Faraday's  law  will  naturally  yield  the  value  of  HZ  at 
this  interface.  As  an  important  example,  for  the  case  where  d  «  0,  Equation 
[16]  will  yield  zero  values  of  both  tangential  components  of  E  and  thus  zero 
values  of  HZ. 

Finally,  the  expressions  for  updating  the  tangential  magnetic  field  compo¬ 
nents  (H  and  H  )  below  the  conductor  must  be  slightly  modified  to  reflect  the 
*  / 

fact  that  the  conductor  sheet  is  protruding  into  the  cells  under  the  slab.  In 

evaluating  Faraday's  law  around  the  contour  (see  Figure  6),  it  is  of  course 

noticed  that  the  point  at  which  Hx  is  to  be  evaluated  is  not  in  the  center  of 

the  contour,  giving  rise  to  an  ambiguity  in  its  value  since  the  slope  of  Hx  at 

that  point  cannot  be  exactly  known.  Fortunately,  since  this  point  is  close  to 

a  perfectly  conducting  surface,  the  slope  of  Hx  with  respect  to  z  will  be 

small  and  thus  relatively  unimportant.  Thus,  the  algorithm  for  calculating  Hx 

needs  only  be  modified  by  noting  that  the  surface  area  of  the  surface  surroun- 

2d 

ded  by  is  smaller  by  the  linear  multiplier  (1  -  —  )  . 
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VI.  CONDUCTOR-BACKED  DIELECTRIC  SLABS  WITH  CRACKS 


The  problem  of  modeling  a  conductor-backed  dielectric  slab  with  a  thin 
crack  in  the  dielectric  presents  a  modeling  challenge  when  using  any  elec¬ 
tromagnetic  technique,  and  the  FDTD  is  no  exception.  This  difficulty  arises 
from  the  relatively  complex  field  distributions  in  the  vicinity  of  the 
crack.  These  fields  exhibit  a  relatively  severe  localized  perturbation  from 
their  relatively  smooth  behavior  away  from  the  crack  due  to  the  collection  of 
bound  charges  across  the  crack.  These  bound-charge  distributions  may  be 
thought  of  as  arising  from  the  disruption  of  the  bound  currents  flowing  within 
the  dielectric,  since  the  crack  presents  an  open  circuit  to  these  currents. 

A  FDTD  model  of  a  conductor-backed  dielectric  slab  with  a  long,  narrow 
crack  is  shown  in  Figure  8.  Here,  it  is  assumed  that  the  air/dielectric  in¬ 
terface  lies  along  a  grid  coordinate  axis,  and  the  conductor  lies  below  at  a 
A  z 

depth  d  (d  <  —j  )  •  The  crack  is  assumed  to  be  a  "clean"  cut  in  the  dielec¬ 
tric,  of  width  w  (w  <  ,  and  extends  through  to  the  conductor.  It  is 

further  assumed  chat  the  crack  runs  out  of  the  page  and  is  aligned  such  that 
it  appears  in  between  grid  lattice  points.  This  alignment  is  such  that  the 
top  of  the  crack  coincides  with  an  evaluation  point  for  E^.  This  particular 


FIGURE  8.  The  orientation  of  a  thin,  conductor-backed  dielectric 
slab  with  a  long,  narrow  crack  in  the  FDTD  grid,  shotring 
the  position  of  the  crack  with  respect  to  the  integration 
contour,  Cy,  used  to  advance  the  tangential  magnetic 
fields  along  the  crack. 
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orientation  of  the  crack  within  the  FDTD  lattice  has  been  chosen  so  that  the 
effect  of  the  crack  can  be  most  accurately  modeled  for  situations  where  the 
incident  field  produces  strong  fields  across  the  crack  (i.e.,  the  dominant 
scattering  mode  of  a  crack). 

Since  the  dielectric  is  assumed  to  be  non-magnetic,  it  can  be  assumed  that 
the  magnetic  field  in  its  vicinity  vary  linearly  throughout  each  cell.  On  the 
other  hand,  this  assumption  cannot  be  made  about  the  tangential  electric  field 
in  the  vicinity  of  the  crack  since  it  would  be  expected  that  this  field 
component  will  be  enhanced  inside  the  crack  and  exhibit  a  distribution  similar 
to  a  static  dipole  in  the  vicinity  of  the  crack,  particularly  within  the  crack 
and  just  above  it.  This  gives  rise  to  two  difficulties.  First,  the  nonlinear 
behavior  will  change  the  way  that  the  right-hand  side  of  Ampere's  law  (Equa¬ 
tion  [la])  is  evaluated  to  calculate  the  field  at  the  lattice  point  given  a 
knowledge  of  its  surface  integral.  A  second  difficulty  stems  from  the  fact 
that  it  is  not  the  value  of  EY(I,J,K)  at  the  lattice  point  that  is  needed  in 
Faraday’s  law  when  updating  magnetic  fields,  but  rather  its  average  over  the 
length  of  its  respective  cell  wall.  Thus,  the  nonlinear  behavior  of  E^  along 
the  £  direction  must  also  be  correctly  modeled  in  order  to  accurately  calcu¬ 
late  the  magnetic  fields. 

The  way  in  which  this  nonlinear  behavior  of  E^  above  and  within  the  crack 
is  modeled  presents  the  most  challenging  aspect  of  this  problem.  As.  is  typi¬ 
cal  in  FDTD  analysis,  there  are  a  number  of  ways  in  which  this  could  be  accom¬ 
plished.  In  the  paragraphs  to  follow,  the  method  that  appears  to  show  the 
most  promise  will  be  presented. 

Since  the  crack  is  assumed  to  be  narrow,  the  capacitance  across  the  gap 
will  be  quite  high.  Thus,  the  polarization  current  flow  in  the  dilectric 

—  + 
(J  *  (e  -e  )  -r— )  will  be  relatively  unaffected  on  either  side  of  the  crack, 
p  o  Ot 

Given  this,  the  crack  can  be  thought  of  as  perturbing  the  fields  of  a  contin¬ 
uous  dielectric  slab  by  subtracting  from  the  total  polarization  current  in  the 


Two  physical  explanations  can  be  offered  to  support  this  claim.  The  first 
is  that  replacing  a  section  of  wire  with  a  large  capacitor  will  have  little 
effect  on  the  current  in  the  rest  of  the  circuit.  The  second  is  that  the 
electric  field  in  the  air  gap  between  the  plates  of  a  capacitor  is  only  mar¬ 
ginally  affected  when  a  thin  dielectric  (much  thinner  than  the  air  gap)  is 
added  between  the  plates. 
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dielectric  that  which  is  no  longer  present  in  the  crack.  This  suggests 
defining  the  total  electric  field  in  the  vicinity  of  the  crack  as: 

E  -  E  -  E  ,  [22] 

a  P 

where  is  the  electric  field  without  the  crack  perturbation  (which  will 
henceforth  be  called  the  unperturbed  field  component),  and  is  the  electric 
field  radiated  by  the  missing  polarization  currents  above  the  ground  plane 
(but  not  in  the  presence  of  the  dielectric).  The  logical  sequence  of  Equation 
[22]  is  depicted  in  Figure  9. 


FIGURE  9.  The  physical  representation  of  Equation  [22],  whereby  the 
fields  due  to  the  long,  narrow  crack  in  a  conductor- 
backed  dielectric  sheet  are  approximated  as  those  with  no 
crack  present,  minus  the  fields  of  the  "missing" 
polarization  currents. 

The  reason  why  the  above  scheme  is  attractive  is  that  the  field  E,  can  be 

d 

calculated  by  the  usual  conductor-backed  dielectric  formulas  (Equation  [17]), 
and  Ep  can  be  calculated  from  the  well-known  field  properties  of  short  ele¬ 
ments  of  current.  Also,  in  the  near  field  of  a  short  current  element  the 
electric  field  is  proportional  to  the  time  integral  of  the  current.^  Thus, 
we  can  write: 


EP  “  /  [Jp  F  +  Jp  F  z<y>]  dt 

v  *y  pz 


where  F  y(y)  and  F  (y)  are  the  near  field  pattern  functions  along  the 


[23; 
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f  direction  of  short  current  dipoles  oriented  in  the  £  and  t  directions, 

respectively.  Jp  and  J  are  the  y  and  z  components  of  the  mis_sing  polari- 

y  2  -  dEd 

zation  currents,  respectively.  However,  since  J  ■  (e  -  £q)  can  be 

written  as: 


E  *  (e  -  e  )  E  F  y(y)  +  (e  -  e  )  E.  F  Z(y) 
p  0  dy  dz 


[24] 


where  E^  and  E^  are  the  and  2-directed  unpreturbed  electric  fields  within 
y  z 

the  crack. 

Knowing  the  relationship  between  the  fields  radiated  by  the  "missing" 
polarization  currents  in  the  dielectric  crack  and  the  relatively  unperturbed 
fields  within  the  dielectric  on  either  side  of  the  crack,  one  can  obtain  the 
following  expression  for  the  ^-directed  electric  field  by  substituting 
Equation  [24]  into  Equation  [22]: 


<Ed  +  A  •  y)  -  Ed  Fy  (y) 


-  Ed  Fy  (y) 
z  J 


[25] 


V  z 

where  F^  (y)  and  F^  (y)  are  the  y  components  of  the  pattern  functions  along 
the  y  coordinate  due  to  and  2-directed  polarization  currents,  respectively, 
and  "A"  is  a  constant.  Notice  in  this  equation  that  the  first  term  is  a 
linear  function  of  y,  but  the  remaining  two  terms  are  nonlinear  dipole  fields. 

Now  that  the  electric  field  along  the  air/dielectric  interface  in  the  vi¬ 
cinity  of  the  crack  has  been  described  in  Equation  [25],  the  line  integral  of 
Faraday's  law  (Equation  [lb])  can  now  be  evaluated.  Looking  at  Equation  [25], 
it  is  obvious  that  the  line  integral  of  the  first  term  will  be  trivial  since 
this  field  component  varies  linearly  along  the  path.  The  integrals  of  the 
remaining  two  terms  will  be  more  challenging. 

One  approximation  of  these  integrals  could  be  obtained  by  lumping  all  of 
the  missing  polarization  currents  at  the  center  of  the  crack.  This  assumption 
would  yield  simple  expressions  for  the  line  integral  of  these  terms.  It  may, 
however,  be  necessary  to  consider  these  currents  as  the  volume  currents  which 
they  truely  are  (i.e.,  they  are  functions  of  depth  into  the  crack).  If  this 
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is  the  case,  the  integration  of  these  terms  will  be  more  complex.  In  either 
case,  however,  it  must  be  remembered  that  the  fields  of  a  short  current  are 
relatively  simple,  so  these  integrals  need  not  be  particularly  complex. 

Summarizing,  this  section  has  proposed  a  simple  model  of  the  behavior  of 
the  tangential  electric  field  in  the  vicinity  of  a  crack  in  a  conductor-backed 
dielectric  slab  that  can  be  calculated  in  terms  of  variables  readily  available 
from  the  FDTD  technique  for  conductor-backed  dielectric  slabs.  From  this 
model  of  the  electric  field  behavior,  a  technique  for  evaluating  the  magnetic 
fields  above  the  air/dielectric  interface  has  been  proposed.  Since  the  dipole 
nature  of  the  fields  in  the  vicinity  of  the  crack  are  built  into  this  tech¬ 
nique,  it  is  expected  to  be  capable  of  accurately  modeling  the  backscatter 
from  such  cracks. 


25 


VII.  CONCLUDING  REMARKS 

This  report  has  concerned  itself  with  the  development  of  a  numerical  tech¬ 
nique  capable  of  calculating  the  backscatter  from  dielectric  slabs,  with  and 
without  cracks.  This  work  represents  an  extension  of  the  well-known  and  pop¬ 
ular  FDTD  technique,  which,  if  not  modified,  is  not  well  suited  to  geometries 
with  fine  detail  (e.g.,  thin  dielectrics  or  cracks). 

In  the  case  of  dielectric  slabs  (with  or  without  conductor  backing),  al¬ 
gorithms  have  been  developed  for  direct  insertion  into  any  standard  FDTD 
code.  These  algorithms  are  similar  in  form  to  the  "standard"  field-advancing 
equations  in  FDTD  codes,  but  allow  for  the  change  in  e  at  the  air/dielectric 
interface  and  the  discontinuity  of  the  normal  electric  field  components  at  the 
interface.  They  are  suited  for  situations  in  which  the  dielectric  slab  is 
small  with  respect  to  wavelength  and  the  cell  dimensions.  Initial  calcula¬ 
tions  from  these  formulas  indicate  that  they  are  indeed  capable  of  advancing 
the  fields  in  the  presence  of  these  scatterers,  although  full-scale  numerical 
testing  will  be  the  subject  of  follow-up  efforts. 

The  analysis  of  conductor-backed  dielectric  structures  with  thin  cracks 
presented  in  this  report  has  identified  the  special  modeling  challenges  of 
this  geometry.  The  formalism  developed  followed  from  a  perturbational 
analysis  of  the  crack  under  the  assumption  that  the  presence  of  a  thin  crack 
has  little  effect  on  the  electric  fields  (thus  the  polarization  currents) 
inside  the  dielectric  on  either  side  of  the  crack.  This  assumption,  along 
with  the  known  radiating  properties  of  small  current  filaments,  yields  a  set 
of  equations  that  depend  on  the  fields  of  the  dielectric  without  the  crack. 
The  best  way  of  evaluating  the  spatial  integrals  in  these  formulas  will  again 
be  the  subject  of  follow-up  work. 
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RADC  plans  and  executes  research,  development,  test  and  selected 
acquisition  programs  in  support  of  Command,  Control,  Communications 
and  Intelligence  ( C 3 I)  activities.  Technical  and  engineering  support  within 
areas  of  competence  is  provided  to  ESD  Program  Offices  (POs)  and  other 
ESD  elements  to  perform  effective  acquisition  of  C3!  systems.  The  areas 
of  technical  competence  include  communications,  command  and  control, 
battle  management,  information  processing  surveillance  sensors, 
intelligence  data  collection  and  handling  solid  state  sciences, 
electromagnetics,  and  propagation,  and  electronic,  maintainability,  and 
compatibility. 


